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FOREWORD 


This  note,  which  presents  results  from  parameteric  flutter 
studies  on  a  typical  section  airfoil  using  piston  theory,  was 
prepared  by  the  Aeroelastic  and  Structures  Research  Laboratory, 
Massachusetts  Institute  of  Technology,  Cambridge  39? Massachusetts, 
for  the  Aircraft  Laboratory,  Wright  Air  Development  Center, 
Wright-Patterson  Air  Force  Base,  Ohio.  The  work  was  performed 
at  the  MIT  under  the  direction  of  Professor  Holt  Ashley  and 
supervised  by  Mr.  G.  Zartarian  and  Mr.  W.  H.  Weatherill.  The 
research  and  development  work  was  accomplished  under  Air  Force 
Contract  No.  AF33 (616) -2482,  Project  No.  1370,  "Aeroelasticity, 
Vibration  and  Noise”,  and  Task  No.  13478,  "Theoretical  Super¬ 
sonic  Flutter  Studies”.  Mr.  Walter  J.  Mykytow  of  the  Dynamics 
Branch,  Aircraft  Laboratory  is  task  engineer.  This  note  con¬ 
tains  data  only  and  represents  a  part  of  a  research  program  on 
the  flutter  of  aircraft  structures  at  supersonic  speeds.  Further 
reports  will  be  published  as  the  research  continues.  The  basic 


research  was  started  1  July  1954  jmd  is,  continuing.  As  part  of 

AD')io  sr*?/- 

this  research  WADC  Technical  Report  56-97,  (Confidential  Report  - 
Unclassified  Title)  "Theoretical  Studies  on  the  Prediction  of 


Unsteady  Supersonic  Airloads  on  Elastic  Wings”  has  been  publish¬ 
ed  in  two  separate  parts.  Part  I,  (Unclassified  ICitle)  "Investi¬ 
gations  on  the  Use  of  Oscillatory  Supersonic  Aerodynamic  Influence 
Coefficients”  which  presents  the  studies  on  the  use  of  the  aero¬ 


dynamic  Influence  coefficient  method,  was  issued  in  December  1955; 

4J i  *  li  ®  v  1 

and  Part  II,  (Unclassified  Title)  "Rules  for  Application  of 

Oscillatory  Supersonic  Aerodynamic  Influence  Coefficients”  which 


presents  working  rules  and  recommendations  for  flutter  analysis 
using  the  aerodynamic -inf luence-coefficient  method,  was  issued 


in  February  1956. 


FOREWORD  (Continued) 


The  authors  are  indebted  to  Professor  Holt  Ashley  and 
Dr.  P.  T.  Hsu  for  their  contributions  to  the  research.  Also 
acknowledgements  are  due  to  Mrs.  Ruth  Lyon  and  her  computing 
group  for  their  help  in  making  the  necessary  calculations, 
Mr.  John  McHugh  for  help  in  preparing  the  tables  and  figures 
and  Miss  Kathryn  Roberts  for  typing  the  final  manuscript. 


ABSTRACT 


An  extensive  set  o£  tables  is  presented  from  binary  flutter 
calculations  on  the  typical  section  airplane  wing  model  of 
Theodorsen  and  Garrick  (Ref.,  1)^  The  airloads  are  predicted 
using  two-dimensional  piston  theory.  The  parameters  studied  in¬ 
clude  thickness  as  well  as  center  of  gravity  position,  elastic 
axis  position,  wing  and  aileron  radii  of  gyration,  and  frequency 
ratio.  The  three  binary  flutter  cases  for  which  data  is  given 
art  bending- tors ion,  bending-aileron,  and  torsion-aileron  flutter. 
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SECTION  1 


PARAMETRIC  STUDY  OF  FLUTTER 
AT  HIGH  SUPERSONIC  SPEEDS 


1.  Introduction 


The  objective  of  this  report  is  to  present  the  numerical 
results  of  extensive  binary  flutter  calculations  on  the  typical 
section  airplane  wing  model  of  Theodorsen  and  Garrick  (Ref.  1) 

'with  two-dimensional  piston  theory  used  for  the  aerodynamic  terms 
in  the  flutter  equations.  The  effects  of  all  dimensionless 
system  parameters  on  the  flutter  eigenvalues  are  included. 

In  order  to  expedite  publication,  no  graphical  presentations 
or  interpretation  of  results  are  given  herein.  This  is  therefore 
an  interim  report,  designed  to  make  the  data  generally  available 
while  a  more  complete  paper  Is  being  prepared  on  the  same  subject. 
The  final  report  will  contain  no  additional  numerical  information 
but  will  be  devoted  to  interpretation  of  findings  and  discussion 
of  ways  in  which  the  material  can  be  used. 

Piston  theory  is  an  approximate  technique  for  estimating 
airloads  on  airfoils  in  unsteady  motion  at  higher  supersonic  and 
low  hypersonic  speeds.  It  had  its  origins  in  work  of  Hayes 
(Ref,  2)  and  Lighthlll  (Ref.  3)  and  has  been  applied  to  aero- 
elastic  problems  in  several  recent  publications  (e.g.,  Refs.  4 
ar.d  5)*  Its  use  for  typical  section  flutter  analyses  is  desirable 
for  several  reasons.  First,  it  yields  a  much  simplified  set  of 
aquations  of  motion,  as  compared  with  more  conventional  unsteady 
aerodynamic  theory.  In  some  cases  the  flutter  determinant  may 
even  be  solved  in  closed  form.  Extensive  flutter  calculations 
using  piston  theory  can  be  made  on  any  small  digital  computing 
machine,  such  as  the  Burroughs  E101  which  was  employed  in  the 

f  resent  Investigation,  Second,  the  influence  of  airfoil  thickness 
and  other  nonlinear  effects,  if  significant)  can  be  Introduced 
without  excessive  additional  computational  difficulty.  To  the 
authors'  knowledge .  this  is  the  first  parametric  study  of  flutter 
which  includes  thickness  effects.  Although  by  no  means  exhaustive, 
the  rentes  of  parameters  covered  appear  to  be  those  of  general  in¬ 
terest  for  the  supersonic  flight  regime. 


Manufcript  released  by  the  authors  May  1957  for  publication  as  a 
Wm  Technical  Mote* 


i  F  * 


In  setting  up  the  flutter  problem  for  this  report,  a  three- 
degree-of -freedom  determinant  was  first  derived  for  the  system 
pictured  in  Fig,  1.  Only  the  three  sub-cases  of  binary  flutter 
were  calculated,  however.  Thus  the  tables  list  data  only  on 
bending- tors ion,  bending-aileron  and  tors ion -aileron  flutter  of 
the  typical  lection. 


The  flutter  equations  and  their  solution  follow  the  pattern 
established  by  Garrick  and  Rubinow  (Ref.  6).  A  symmetrical 
double-wedge  airfoil  shape  has  been  chosen  as  typical  of  super¬ 
sonic  configurations. 


Section  2  gives  the  computational  details,  while  Section  3 
discusses  the  manner  of  tabular  presentation  of  data.  In  this 
latter  connection,  a  few  remarks  are  in  order.  Past  convention 
has  usually  been  to  present  compressible-flow-flutter  boundaries 
on  plots  of  dimensionless  speed'  (e.g.  )  vs.  Mach  number  M. 

On  such  a  plot,  flight  of  a  given  aircraft  at  a  given  altitude  is 
represented  by  a  sloping  straight  line  through  the  origin  of 
coordinates.  A  flutter-free  condition  for  the  given  altitude 
occurs  at  any  M  where  the  actual  flutter  curve  lies  above  this 
"altitude  line".  This  form  of  presentation  has  the  disadvantage 
that  the  airspeed  is  contained  implicitly  in  both  the  ordinate  and 
abscissa  of  the  stability  diagram.  It  has  been  pointed  out 
(Ref.  6)  that  the  flutter  stability  boundary  can  equal!"  well  be 
presented  as  a  curve  of  fa<u%,vs.  M,  where  a  is  the  spec  t  sound 
at  flight  altitude.  Flight  of  a  given  aircraft  at  a  given 
altitude  is  then  described  by  a  horizontal  straight  line;  no 
flutter  occurs  if  this  straight  line  lies  above  the  actual  flutter 
curve.  Alternatively,  the  flutter  curve  can  be  interpreted  in 
terms  of  the  torsional  stiffness  required  for  safe  flight  of  a 
wing  of  given  geometry  at  a  given  speed  of  sound  and  Mach  number. 

It  has  been  pointed  out  that  the  curves  for  various 
altitudes  can  be  brought  in  closer  juxtaposition  by  multiplying 
the  ordinate  bu>-A  by  some  power  of  the  wing-to-air  density  ratio 
(thus,  can  be  interpreted  in  terms  of  indicated  rather 

than  true  airspeed).  For  this  reason,  and  because  and  M 
always  appear  as  the  product  in  the  analysis,  all  flutter 
data  in  the  present  report  are  given  as  tables  of  /■*•%  (orSdt%  ) 
vs.  /a.M,  each  row  representing  the  ordinates  of  the  stability 
boundary  for  one  set  of  dimensionless  parameters. 


*See  List  of  Symbols  for  definitions  of  all  physical  quantities  * 
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2.  Development  of  Flutter  Equations 


Figure  la  illustrates  the  two-dimensional  dynamic  model 
which  forms  the  basis  of  this  trend  study  on  supersonic  flutter. 

It  is  identical  with  the  ‘‘typical  section”  model  used  by 
Theodorsen  and  Garrick  (Refs.  1  and  6).  This  system  has  three 
degrees  of  freedom,  described  by  the  time -dependent  quantities 
h,  ot  ,  and  £  ,  which  are  analogous  to  bending,  torsion  and  con¬ 
trol-surface  rotation,  respectively,  of  an  unswept  wing  or  tail 
surface.  When  the  equations  of  motion  of  this  model  in  a  super¬ 
sonic  airstream  are  set  up,  under  the  assumption  that  h,  ot  ,  £  , 
vary  simple -harmonically  with  time,  there  results  the  well-known 
flutter  determinant  (cf.  Ref.  6).  This  complex  determinant,  whose 
vanishing  supplies  two  characteristic  equations  for  the  flutter 
eigenvalues  UF  and  4>,  ,  is  of  third  order  in  the  present  example: 

ABC 

D  E  F  *  0 

G  H  1  (1) 

In  terms  of  quantities  defined  in  the  List  of  Symbols,  the 
elements  of  Eq.  (1)  are  the  following: 

B  =  -/JLX«  +  L3  +  iLj 

c-  -fXf-fL'-nLi 

D  «  ~ /A** 

E  -  /trffo-l} 

F  "  [rf  +2(*.  -x.)Xj7 

G  * 

H  -  -yu/V  +2(x, -*•)*,!  +  #1+1*4 

1  *  /-  il ~'l  *  ** tW* 

(2) 

Here  (  )  >  *tc.,  are  dimensionless  representations  of 

the  lifts,  pitching  moments  rat  d  control -surface  hinge  moments  due 
to  simple  harmonic  motions  in  the  three  degrees  of  freedom.  The 
definitions  adopted  here  coincide  with  those  originally  given  by 
Garrick  and  Rubinow  (Ref.  6) .  For  sufficiently  high  supersonic 


WADC  TH  57-310. 


3 


- - 2  bx, - • 
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(o)  DYNAMIC  MODEL  OF  AIRFOIL  SECTION 


(b)  GEOMETRICAL  SHAPE  OF  AIRFOIL  SECTION 

Fig.  1  The  Model  Used  For  Flutter  Analysis 
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flight  Mach  number  M,  these  coeffients  can  be  put  into  the 
relatively  simple  forms  yielded  by  piston  theory  (Ashley  and 
Zartarian,  Ref.  5).  Assuming  a  syoaetrical  double-wedge  airfoil 
section  of  chord  2b  and  thickness  ratio  &  ,  and  accounting  for 
the  second-order  effect  of  thickness,  Zartarian,  Heller  and 
Ashley  (Ref.  4)  give  the  following: 

L,  =  0 


H  =  7T* /('-■** J 

l6  -  *&{•-*/**} 


*2 

M4  -  jfc  {j 

I /  -  2^,Ms] 


0-xM-2*.)+xJ  ! 

HA*  V 

(l  -H,)Xf/l~2*m)+$‘*i 

HA 

0 

l!=g{l-2A*l] 


Wt‘ 

40-k y. 


’-*/**} 


(3) 


As  in  Ref.  (1),  the  present  report  concentrates  on  the 
influence  of  various  system  parameters  on  solutions  of  the 
flutter  determinant  MEq  (l)“  for  the  sub-cases  obtained  by 
successively  suppressing  each  of  the  degrees  of  freedom.  The 
determinants  for  these  three  binary  flutters  are 
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(1)  Bending  -  torsion  flutter: 

A  B 

=  0 

D  E  (4) 

(2)  Bending -aileron  flutter: 

A  C 

-  0 

G  I  (5) 

(3)  Tors ion -aileron  flutter: 

E  F 

-  0 

HI  (6) 

It  is  significant  that  the  bending- torsion  sub -case  can  be 

solved  explicit^  for  the  dimensionless  speed  and  frequency  of 
flutter: 


These  expression*  display  explicity  the  influence  of  the  various 
system  parameters:  and  combinations  thereof  on  the  eigenvalues. 

It  is  unfortunate  that  similar  closed-form  solutions  cannot  be 
mritten  down  for  the  other  tvo  binary  flutter  cases.  From  the 
computational  standpoint,  however,  the  implicit  equations  yielded 
by  the  bending-aileron  and  tors  ion -aileron  determinants  can  be 
solved  for  speed  and  frequency  with  nearly  equal  ease. 

3*  Description  of  Results 

The  combinations  of  parmieters  used  for  each  binary 
flutter  solution  are  listed  in  Tables  1-3. 

Tables  4  through  62  contain  all  the  calculated  results. 
Hank  spaces  indicate  that  nc  calculations  were  mede.  Asterisks 
indicate  that  careful  examination  of  the  equations  shows  that  no 
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solutions  exist.  Dashes  in  the  boxes  mean  that  calculations  were 
made  and  no  real  solutions  were  found. 

The  numbers  contained  in  the  tables,  as  mentioned  in  the 
Introduction,  are  values  of  nAd^A  at  flutter  for  the  particular 
combination  of  parameters  involved.  The  tables  are  setup  with 
the  purpose  of  making  curves  of  yuht^/a.  versus  yuM  for  a  series  of 
values  of  thickness  ratio  with  other  parameters  fixed.  was 

chosen  because,  for  a  given  configuration  (values  of  h,  a)*  and  m 
fixed)  and  altitude,  is  represented  by  a  straight  hori¬ 
zontal  line  in  a  graph  of  versus  ^w-M.  Also,  a  flight 

diagram  showing  the  design  capabilities  of  a  particular  aircraft 
may  be  superimposed  on  a  plot  of  versus  yC<.H.  The  Mach 

number,  M,  always  appears  in  linear  combination  with  yu  ,  thus 
was  considered  to  be  the  most  convenient  and  economical  scale  for 
the  abscissa  (ef.  Eq.  7)*  Note  that  ya  also  appears  with  the 
thickness  parameter  in  the  single  comoination  .  An  example  of 
drawing  a  flutter  boundary  on  a  flight  diagram  is  given  in  the 
Appendix.  Although  the  tables  are  set  up  to  oe  used  to  show 
trends  due  to  thickness,  the  reader  should  have  no  trouble  in 
combining  any  set  of  numbers  to  illustrate  desired  trends. 

All  the  results  were  obtained  for  a  typical  section  of  a 
symmetrical  double  wedge  airfoil  section  of  chord  2b  and  thickness 
2&b,  For  the  modes  used  in  this  report  (first  power  in  x  only), 
as  shown  in  ref.  5,  the  actual  profile  of  the  wing  enters  the 
piston  theory  aerodynamic  derivatives  only  in  the  form  of  the  area 
of  the  airfoil  section,  the  area  of  the  flap,  the  first  area  mo¬ 
ments  of  the  wing  and  flap  about  their  leading  edges,  and  the 
thickness  of  the  profile  at  the  hinge  line,  all  non- dimens ionalized 
with  respect  to  the  semichord.  If  these  parameters  can  be  made 
numerically  equal  for  two  different  profiles  (say  by  adjusting  the 
thickness  parameter),  then  the  two  profiles  would  have  the  same 
aerodynamic  derivatives.  This  may  be  done  only  for  special  systems 
for  the  binary  cases  including  the  aileron  mode  and  for  the  ternary 
system  where  three  or  more  of  the  above  conditions  must  be  satisfied. 
However,  for  the  bending- tors  ion  system,  where  only  the  area  and 
first  area-moment  conditions  must  be  satisfied,  similarity  rules  may 
be  found  between  all  profiles  having  the  same  non-dimensional  cen¬ 
troid  of  area  location.  For  example,  a  doubly  symmetric  biconvex 
airfoil  section  with  parabolic  thickness  distribution  of  maximum 
thickness  2b  6^  will  have  the  same  bending -tors ion  flutter  bound¬ 
ary  as  the  double-wedge  section  of  this  report  of  maximum  thickness 
2b  (  ‘t&Mc/s ) .  This  of  course  assumes  that  all  other  system  para¬ 
meters  ^re  identical.  These  various  aspects  will  be  further  discuss¬ 
ed  in  the  final  report. 
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No  real  solutions  £0;  values  of  Ok.  Above  0.8  for 
the  range  of  fi m  considered.  ^ 
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TABLE  56  £ _ 1  FOR  BENDING-AILERON  FLUTTER 
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APPENDIX 


The  use  of  the  parameters  and  /tM  to  plot  a  flutter 

boundary  enable  one  to  plot  the  requirements  of  the  aircraft  and 
the  flutter  boundary  on  the  same  graph.  For  example,  Fig.  2 
shows  the  region 


Fig.  2 

Graph  of  Flight  Capability  of  Typical  Aircraft 
in  Terms  of  Altitude  and  Mach  Number 


of  expected  flight  capability  for  a  typical  aircraft  on  an 
altitude  versus  Mach  number  plot.  After  setting  the  wing  config¬ 
uration  (which  includes  4  ujk  and  m0)  of  the  aircraft  involved. 
Fig.  2  may  be  redrawn  on  a  versus  yuM  plot  and  should 

look  something  like  the  sketch  in  Fig.  3. 


Fig.  3 


Flight  Capability  of  Typical  Aircraft 
in  Terms  of  and  yuk^^a 
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Finally*  if  we  superimpose  the  flutter  boundary  on  the  graph  of 
Fig.  3 »  we  will  have  a  sketch  similar  to  the  one  shown  in  Fig.  4. 
This  diagram  immediately 


STABLE 


Fig.  4 

Flutter  Boundary  Superimposed  on  Flight  Boundary 
for  Typical  Aircraft  in  Terms  of  uM  and 


shows  what  part  of  the  required  flight  region  of  the  specific 
airplane  configuration  is  unstable. 
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